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Abstract 

We present model-independent measurements of the branching fractions B{B^ D*^D*^) and 
B{D+ (/>7r+) based on 19.3 fb~^ of data collected by the BABAR detector at the PEP-II 
e~^e~ B Factory. Neutral i?-meson decays to the D*~^D*~ final state are selected with a partial 
reconstruction of the D*^; that is, only the D*~ and the soft photon from the decay D*^ — > Df'y 
are reconstructed. The branching fraction B{B^ D*^D*~) is extracted from these event yields, 
while B{Df — > (j)TT~^) is determined by combining this result with a previous measurement of the 
product B{B^ D*~^D*~) x (/)7r+) with partial reconstruction of the D*~ . We obtain 

the following preliminary results: 

B{B^ D*+D*~) = (1.50 ± 0.16 ± 0.12)%, 
B{Dt (j)Tr+) = (4.7 ± 0.6 ± 0.8)% 

where the first error is statistical, and the second systematic. 
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1 INTRODUCTION 



We present a measurement of the branching fractions B{B^ D*^D*~) and 13{Df — > cp-K^) using 
a partial-reconstruction technique^ A precise measurement of the branching fraction for this 
mode is important because nearly all Df branching fractions are determined by normalizing the 
measurements to B{Df — > (j)TT~^) [2]. The present uncertainty on B{Df — > (f)TT^) thus affects many 
of the results regarding mesons, including the determination of the decay constant by means 
of purely leptonic decays and the measurement of the D'^ — > K inclusive decay rate, as well as 
6-physics analyses where a Df or a D^'^ is fully reconstructed. 

In the factorization model for two-body decay rates, it is assumed that each contribution to the 
transition amplitude of the process is the product of two currents that can be evaluated separately. 
This model has been successful in describing the measured branching fractions and polarizations 
for B meson decays such as 5° D*-tt+ X) ^ D*' p+ and D*-af 0, in which the 

momentum transfer in the process is low (g^ ~ M^,Mp). Measurements of decay rates for modes 
such as B^ D*^D*~ (Fig. EJa)) allow tests of the predictions made using the factorization 
model when the W emits a light and a heavy quark and so the momentum transfer is high {q^ ~ 

The Feynman diagram for the decay B{Df </)7r+) is shown in Fig. ^b). 




2 THE BABAR DETECTOR AND DATASETS 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage ring 
and correspond to an integrated luminosity of 19.3 fb^^. A detailed description of the detector can 
be found in Ref. jTj. 

In addition to this data sample, several simulated event samples were used in order to study 
efficiency and backgrounds. For background studies, we used Monte Carlo samples of B^B^ events 
(equivalent to an integrated luminosity of 270 fb~^), B^B^ (150 fb^^), e+e"^ cc (70fb~^) and 
e+e~— > uu, dd, ss (70fb~^). We used two signal samples in which the B'^ — > D*'^D*~ decay 
proceeds either with completely longitudinal or transverse polarization; an additional signal sample 

^Here and in the following, charge-conjugate processes are implicitly considered. 
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was extracted from the B^B^ by selecting only B^ D^^D* decays, with no further restriction 
on the D*~ and decays. 

3 ANALYSIS METHOD 

3.1 Analysis Strategy 

The partial reconstruction technique results in a significantly larger sample of events than full 
reconstruction. Moreover, the B^ — > D*~^D*~ — > {Dfj){D^TT~) decay is interesting from an 
experimental point of view since it can be partially reconstructed in two ways: the D*~^ can be 
fully reconstructed and combined with the slow pion from the decay D* D^tt^ , or the D*~ can 
be fully reconstructed and combined with the soft photon from the decay D*^ Df^. 

The former technique has been used in BABAR |Sj to measure B{B^ — > D*^D*~). However, the 
precision one can achieve through this technique is limited by uncertainty on B{Df — > (pir'^). By 
applying this method, the B^ — > D*^D*^ branching fraction can be expressed as 

13{B^ ^ D*+D*-) = — ^ ^Dt+n- 



where Nj^*+^_ is the number of partially reconstructed D* candidates, is the number of 

neutral B meson pairs, B^ " are the branching fractions, Si are the total reconstruction effi- 
ciencies^ computed separately for each Df decay mode, and the index i runs over all Df decay 
modes used in the reconstruction (in Ref. [S], Df (/)7r+, Df K*^K°, and Df K*^K^). 
Partial Dl'^ reconstruction similarly yields 

^f^o ^ ^*+^*-^ \ Np,-^ ^ 

^ ' ' 2N^^ BiDt^ ^ Dt-i)B{p*- ^ • ) ' 

where Nj^*-^^ is the number of partially reconstructed D\'^ candidates. The result now depends 
on ■Sj^", the branching fractions for the modes, which are measured much more precisely than 
those of the . 

The — > ^'k'^ branching fraction can be extracted by combining the two methods. Dividing 
Eq.^by Eq.l^Jand solving for B{Df — > 07r+), this last quantity can be expressed as 

where R^" = B^'' /B{Df (f>7r^) is the branching fraction of each Df mode relative to the Df —> 
(f)7r~^ mode, and the D*^ir^ {D*^^) yields are normalized to the same luminosity. In the systematic 
uncertainty determination the contributions given by B{D*~ — > D^7r~) and B{D*~^ — > D^j) are 
clearly cancelled according to Eq. 01 

3.2 Signal Extraction 

We reconstruct the B^ D*^D*^ (Df^){D^Tr~) decay by combining photons in the event with 
fully reconstructed D*^ mesons, without requiring reconstruction of the Df from the D*^ decay. 

^Both exclusive D*'^ reconstruction and partial reconstruction efficiencies are included in e^. 



9 



In order to extract the signal, we compute the missing mass Mmiss recoihng against the D* -7 
system 



Mmiss = Y (-Ebeam " Ed* - E^f - {l^ B - 'f D' " 'f -yY ■ (4) 

For signal events, this must be the mass within experimental resolution. The kinematics of 
the event are not fully constrained with this procedure and one of the decay parameters must be 
chosen in an arbitrary way. In particular, approximating the energy of the B meson in the e+e~ 
center-of-mass (CM) to the CM beam energy, the opening angle between the B momentum vector 
and the measured D* momentum vector can be calculated from 4-momentum conservation in the 
S° ^ D*+D*- decay 

m| - m^,+ - 2EbEd* 

cos 'dBD- = oi-" I • 

^\Pb\\pd*\ 

The B four-momentum is therefore determined up to the azimuthal angle around the D* direction. 
However, an arbitrary choice of the azimuthal angle {cos(pBD* = 0) introduces only a negligible 
spread (less than 1.5 MeV/c^) in the missing mass distribution. 



3.3 Event Selection 

To reject events from continuum, we require the ratio of the second to the zeroth Fox- Wolfram 
moment {R2) (9j to be less than 0.3. 

Candidates for D*^ are reconstructed in the D^ir^ mode, using decays to K^tt^ , i^^vr^vr+vr" 
K~^tt~tt'', and Kg7r~^-K~ , here listed in order of decreasing purity. The probabilities of both the 
and D* vertex fits are required to be greater than 1%. The D* momentum in the T{4S) 
frame must satisfy 1.4GeV/c < p^^^{D*-) < 1.8GeV/c. Moreover, we require the reconstructed 
mass of the particle to be within 2.5 standard deviations of the nominal mass, and the 
D*- Q-value (Q(D*-) = M{D*-) - M{D^) - M(7r-)) to satisfy Qio < Q{D*-) < Qhi, where 
Qio = 4.00 to 5.25 MeV/c^ and Qhi = 6.75 to 8.00 MeV/c^, depending on the D° decay mode. 
Kaon identification is required for the modes K~^7r~7r° and K~^7r~Tr~^7r^ . For the mode KgTr~^7r~, 
the invariant mass of the vr+vr^ from the Kg decay is required to lie within 15MeV/c^ of the Kg 
nominal mass and its flight length must be greater than 3 mm. If more than one D* candidate 
is found, for each decay mode we first select the candidates in which the pion from the decay 
D^TT~ has hits in the drift chamber. Among these, the one with the minimum value of 

= [{Q{D*-) - QpDG{D*-))/^QiD'-)? + m{D^) - MpDG(I)°))Mf(DO)]' is retained. Finally, 
if candidates from different decay modes are present, we select the one with the best average 
purity. 

The selection of photon candidates is based on the optimization of the statistical significance 
of the observed signal {S/\/S + B, where S and B are the number of signal and background 
photons), using generic Monte Carlo events. We apply a vr'^ veto on photon candidates, rejecting 
them if their invariant mass, calculated with any other photon candidate in the event, is between 
115 and 155MeV/c^. The following additional cuts are applied on the photon energy in the T{AS) 
CMS {E'^^^), the cluster lateral moment (LAT) and Zernike moments [n] of order {2, 0} (Z20) 
and {4,2} (Z42): S^^^ > 142 MeV, 0.016 < LAT < 0.509, Z20 > 0.85, Z42 < 0.14. If more than 
one photon is found in the event, we choose the one which maximizes the value of a likelihood 
ratio based on four photon variables {E, E^^^, Ncry, LAT), where E is the photon energy in the 
laboratory frame and Ncry is the number of calorimeter crystals involved in the electromagnetic 
shower. 
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3.4 Selection efficiency and Monte Carlo validation 



To effect tlie partial D^^ reconstruction in Monte Carlo events, the Monte Carlo sample is split in 
two parts. The signal reconstruction efficiency is determined from — > D*'^D*~ events extracted 
from the first half of the sample by performing a minimum-x^ fit to the missing mass distribution. 
The signal peak, centered on the nominal Df mass, is well described by a Gaussian probability 
density function (p.d.f.), while the background, which is mainly due to random D*—j combinations, 

is parametrized with the function B(m) = a (l — e~''('"~™'™=''=)) ( ) , where m = M^iss and 

iTT'ma.x is the end point of the missing mass distribution. We perform a single fit to all decay 
modes; the sum of the branching fraction-weighted efficiencies for the four reconstruction modes 
is computed from the number of signal events fitted in the range \m — M{Df)\ < 41MeV/c^, and 
found to be {eB) = J^j i^j ■ J3f ) = (7.14 ± 0.16) x IQ-^. 

We have validated the fitting technique and the method of extracting the signal on the other 
half of our Monte Carlo sample. The distribution of the missing mass of partially reconstructed 

candidates is shown in Fig. [^Jfor B^B^ (including signal), B^B^ , and continuum Monte Carlo 
events. The signal yield is extracted from a minimum-x'^ fit of the missing mass distribution to a 
sum of the signal, described by a Gaussian function, and the background, described by the p.d.f. 
introduced above. From the signal yield, using Eq.|2l we obtain ^(5*^ D*^D*~) = (1.43ib0.04)%, 
which is consistent with the value of 1.41% used in the generation of the Monte Carlo. 




Figure 2: Monte Carlo simulation of missing mass distributions. The missing mass is defined by 
Eq. 01 Contributions of continuum (black), B~^B~ (grey) and B^B^ (points) are added. The solid 
line shows the fit to the signal plus the sum of all backgrounds. The dashed line is the fit to the 
background component only. 



3.5 Results on data 

Figure shows the missing mass distribution in our data sample. The same fitting procedure 
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Figure 3: Missing mass distribution in the data sample. The sohd hne shows the fit to the signal 
plus background, the dashed line is the fit to the background component only. 



applied in the previous section to the Monte Carlo sample is used to extract the number of signal 
events in the data sample. In the fit we let all parameters float except the mean and the stantard 
deviation of the Gaussian signal, which are fixed to their Monte Carlo values. The result of the fit 
to the missing mass distribution is shown in Fig. |31 as well. The signal yield in the data sample is 
Nr^+ = 1382 lb 145 events. The of the fit is 53.3 for 54 degrees of freedom, corresponding to a 
probability of 50.1%. 

From this yield we obtain B{B^ D*+D*-) = (1.50±0.16)%, where the error is just statistical. 

4 SYSTEMATIC STUDIES 

The main sources of systematic uncertainties on the D*^D*~ branching fraction measure- 

ment are listed in Table ^ The Monte-Carlo-statistics uncertainty is due to the statistical error 
on the efficiency determination. The uncertainty on the background shape is evaluated by fitting 
the missing mass distribution using a different p.d.f. for the background^, and assigning the rel- 
ative signal yield difference as systematics. The systematic uncertainty due to tracking efficiency 
is evaluated applying a random inefficiency of 0.8% per track (1.6% for the soft pions from D* 
decays). We assign as an uncertainty the difference between the yield obtained in this way and 
the one described in Sec. 13. 41 The systematics associated to photon reconstruction efficiency and 
particle identification are evaluated in a similar way. We find a 7% difference in the overall se- 
lection efficiency between our samples with complete longitudinal or transverse polarization in the 
D*^D*~ decay. The uncertainty due to the dependence on polarization is computed taking 
into account the measured value |Hj of the fraction of longitudinal polarization and its uncertainty 

^The alternative baekground p.d.f. has the following functional form: B{m) — — mmax) ^ 

c + (m - mmax)" 
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Table 1: Fractional systematic uncertainties (%) for the D*^D* branching fraction mea- 

surement. 



Source 


Error (%) 


Monte Carlo statistics 


2.3 


Background shape 


2.9 


B counting 


1.1 


Tracking efficiency 


2.4 


Soft pion efficiency 


1.6 


Photon efficiency 


4.2 


Particle identification 


1.5 


Polarization uncertainty 


0.8 


branching fractions 


3.2 


B{D*- Wtt~) 


0.7 


B{D*+ ^ Dti) 


2.7 


Total systematic error 


7.9 



Vl/T = (51.9 lb 5.7)%. Finally, the uncertainties on D^, D* and D^^ branching fractions are 
propagated through the analysis. 

5 PRELIMINARY BRANCHING FRACTION RESULTS 

The B^ D*^D*^ branching fraction is found to be: 

B{B^ D*+D*~) = (1.50 ± 0.16 ± 0.12)%, (6) 

where the first error is statistical, and the second systematic. The Df (pn^ branching fraction 
can be extracted by comparing this result with the measurement of the B^ D*^D*~ decay with 
partial D*' reconstruction ^(5° ^ D*+D*') = (1.97 ± 0.15stat ± 0.30syst ± 0.49^+^^^+)%. 
In this measurement the world average branching fraction B{Df (j)TT^) = (3.6 it 0.9)% was used. 
From Eq.Elwe obtain therefore: 

B{Dt 07r+) = (4.7 ± 0.6 ± 0.8)%, (7) 

where the first error is statistical, the second systematic. The systematic uncertainty on this 
branching fraction is dominated by the measurement using partial D*~ reconstruction. 

6 SUMMARY 

A measurement of the B^ D^^D*^ branching fraction is performed, using data corresponding 
to an integrated luminosity of 19.3 fb^^, with a partial reconstruction technique of B meson: 

B{B^ D*+D*-) = (1.50 ± 0.16 ± 0.12)%. 

This result is compatible with, and improves on the precision of previously published experimen- 
tal results [21 m, and should be compared with the most recent theoretical results based on the 
factorization assumption ^: B{B" D*+D*-)theor = (2.4 ±0.7)%. 
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The (pTT^ branching fraction result is 

l3{Df (Ptt+) = (4.7 ± 0.6 ± 0.8)%. 

This new determination is compatible with the published CLEO result T and a preliminary mea- 
surement from Belle |12j . 
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